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Abstract 

We calculate the two time current correlation function using the experimental data of the current- 
time characteristics of the Gas-DNA-decorated carbon nanotube field effect transistor. The pattern 
of the correlation function is a measure of the sensitivity and selectivity of the sensors and suggest 
that these gas flow sensors may also be used as DNA sequence detectors. The system is modelled 
by a one-dimensional tight-binding Hamiltonian and we present analytical calculations of quantum 
electronic transport for the system using the time-dependent nonequilibrium Green's function 
formalism and the adiabatic expansion. The zeroth and first order contribution to the current 
J*-*^^ (t) and J^-^'^ (t) is calculated where J^^^ (t) is the Landauer formula. The formula for the time- 
dependent current is then used to compare the theoretical results with the experiment. 
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Study of electronic transport through mesoscopic or nanoscale systems is one of the 
most active areas of today's research. Electronic transport in these systems is divided into 
stationary and time-dependent phenomena. The stationary transport for nonequilibrium 
systems has been studied by many authors Meir and Wingreen reformulated the 
ideas given in [l| (and references therein) using the Keldysh approach to study interacting 
mesoscopic systems which led to a Landauer formula. Later, Meir, Wingreen and Jauho 
introduced a general formulation for time-dependent transport through mesoscopic (nano) 
structures j^. Using this formalism theoretical understanding of experiments js], O] has been 
presented for time-dependent voltages. 

In this work, we derive an expression for the current for gas flow over a single-stranded 
DNA (ss-DNA)-decorated single wall carbon nanotube (SWCN) connected to the source (S) 
and drain (D) contacts maintained at zero gate voltage Vg = OV and a fixed bias voltage 
Vft = lOOmV. For this system, time-dependence arises as different number of gas molecules 
flow over and interact with the DNA-decorated SWCN at each time i s for an exposure time 
of 50 s. This results in the hopping integral and on-site energy to be functions of time which 
leads to a time-dependent Hamiltonian for the SWCN. Thus, to deal with gas flow over nano- 
structures the time-dependent nonequilibrium Green's function (NEGF) formalism is well 
suited. The choice of the model is motivated by the experiment |5| which finds the current- 
time characteristics of the ss-DNA-decorated p-type SWCN-FET sensors upon odour and 
air exposures. The odours are dimethyl methylphosphonate (DMMP)^, 2, 6 dinitrotoluene 
(DNT), trimethylamine (TMA), propionic acid (PA), and methanol [5;]. The characteristic 
curves show fluctuations in the current response. We calculate these fluctuations in terms 
of the two time current correlation function. 



I. TWO TIME CURRENT CORRELATION FUNCTION 



We begin with 
experimental data 



;he calculation of the two time current correlation function using the 
sl. The two time current correlation function C{t) is defined as 
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c{i) = -J2I^iw+st) --j2m X -j2h{t+st}, (1) 
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where /j(t) = liit)/!^ is the normalized source-drain current at time t in the i*^ gas exposure 
cycle. Let us consider 5t = 1 s and n is the total number of exposure cycles. 



The correlation function calculated here is different from that shown in our previous 
work equation (2) and for other mesoscopic systems equation (5.57) and js]. The 
second term in equation (2) of [6'] and equation (5.57) of j3] is < Ii{t) >^ and < g >^ 
respectively, where < ■ ■ ■ > denotes the ensemble average. Instead of < /j(t) here 
we use ^"=1 -^j(^) ^ Sr=i + average ^ X^ILi -^^(^ different from the 

average ^ Y17=i hit+^'i)) because for each 5i there is contribution to the current from another 
base. In the system (Gas+DNA+SWCN complex), at 5i = ii s the gas molecule interacts 
with a base (say, cytosine) attached to the carbon atom and for the next 5i the gas molecule 
interacts with another base thymine (say), so for each 5i the current changes as the Gas- 
DNA-base complex changes upon gas exposure. This is not the case for the usual disordered 
mesoscopic systems considered in experiments, where both the averages ^^"=i-^j(^ and 
^ Sr=i + same ( 7| equation(5.57) and [sl). 

We observe distinct patterns for the different odours and sequences, figure 1. This indi- 
cates selective recognition of each odour by the sensors. In particular, the pattern of the 
correlation function for methanol with sequence 1, figure 1(c), is different from methanol 
with sequence 2, figure 1(d). This shows that the correlation function is sensitive to different 
gases and DNA sequences. Hence, the two time current correlation function is a measure 
of the sensitivity and selectivity of the DNA-decorated SWCN sensors. The correlation 
function suggests that these gas flow sensors may also be used as DNA sequence detectors 
where the pattern of correlation functions may be used as a benchmark for the particular 
chemical signal encoded in a DNA sequence. Figure 1 shows highly fluctuating data, these 
fluctuations are due to the structure of the DNA sequences as shown in 



II. TIGHT-BINDING MODEL AND NONEQUILIBRIUM GREEN'S FUNCTION 

The experiment [sl with gas exposure indicates that the sensor response (AI/Iq) is zero 
for pristine SWCN and as DNA is applied on SWCN the sensor response changes. Since 
the ss-DNA sequence is a linear chain we can model the ss-DNA-decorated SWCN using 
a one-dimensional tight-binding Hamiltonian where the electron hops between two carbon 
atoms attached to the DNA bases with different hopping integrals and on-site energies at 
different times. In our formalism, the effect of the gas flow and DNA functionalization in 
the channel can be modelled by a time- dependent potential. Here, we model the effect by 
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the time-dependent hopping integrals and on-site energy in a self-consistence manner {q]. 

Now we investigate the electronic transport through the model system using the time- 
dependent NEGF formalism. For this, a tight-binding model of the SWCN and Gas-DNA 
complex is set up. In this microscopic tight-binding model, there are two approximations: 
the first approximation is the tight-binding between the carbon atoms of SWCN and the 
second is the interaction between the carbon atoms of SWCN and Gas-DNA-base complex 
0, [3], where we assume that the carbon atom interacts with the nearest-neighbor DNA 
base. 

Figure 2 shows the schematic diagram of the chosen model system, i.e., Gas-DNA-SWCN 
sandwiched between S and D contacts. Here, we present a simplified picture of the complex 
model and its operation. Ai, Bi, A2, B2, ■ ■ ■ , Xi represent the carbon atoms, where X can 
be A or B and the index i=l, 2, ■ ■ ■ , M is the number of the carbon atom. The first and 
the last carbon atoms, A-i and Xm are connected to S and D. The bases cytosine, thymine 
e.c. of a„ ..DNA _ 2 ^ (.how„ by ova.) a«ac.ed .o d.ffe.e. ca.bo„ a.on. 
and the circle represents the gas molecule. The arrow indicates the path of transmission. 

The operational principle of the model is based on the changes in the electrical properties 
due to DNA bases and gas molecules adsorbed on the SWCN surface. Initially, the DNA 
sequence is applied on the SWCN, e.g., the cytosine base attaches to the carbon atom Bi 
and thymine base attaches to A2 and so on. When the gas is exposed to the sensor for 
duration of 50 s then at time ii s the gas molecule interacts with cytosine base through vdW 

nri 

forces and/or mutual polarization [6|, |10|. As a result, the electron hops from the vr orbital 
of one carbon atom Ai to the neighboring carbon atom Bi and the tunnelling of the electron 
from Ai to Bi is an elastic process with the corresponding integral referred to as the hopping 
integral 711 (ti). Then, the electron hops from Bi to A2 with the hopping integral 712(^1), and 
from A2 to B2 with 722(^1) = 70? where 70 is the hopping integral without the gas and DNA. 
Hence, the hopping integrals between other carbon atoms 7ij(ti) are also 70, figure 2(a). 
The time scale of electron transport through the SWCN is far less than the experimental 
time scale of the gas flow, therefore the electron sees a steady state at ti s. Similarly, at 
^2 s another gas molecule interacts with the thymine base attached to the carbon atom A2, 
resulting in the hopping integrals 711(^2), 712(^^2), 722(^2), and other 7ij(t2)=7o, figure 2(b). 
In this case, the electron sees a steady state at ^2 with 71^(^2) different from 7ij(ti) as the 
bases interacting with the gas molecules are different. In a similar way at time tM, the time 
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when the gas molecule interacts with the last base of the DNA sequence attached to the 
carbon atom, the hopping integrals are 711 (tAf), 7i2(^a/), ■ ■ ■, Ipqiiu), figure 2(c). Hence, 
we find the hopping integrals as well as the on-site energy change with the time i at which 
the gas molecules trigger the different bases of the DNA sequence attached to the SWCN 
and depend on the DNA bases, gases, and geometry of attachment of the bases and gases to 
the SWCN. Because of the time-dependent hopping integral and on-site energy the relevant 
Green's functions for the SWCN (Gc) are also time- dependent. 

As the experimental time scale of the gas fiow is larger than the time of electron transport 
inside the SWCN (~ lO^^^s) one can express the leading transport properties of the metal- 
SWCN-metal system in terms of the Green's function given by G{t,t) = [t- H{t)]-, where 
e = e ±11] with irj an infinitesimal imaginary term and H{t) is the Hamiltonian at time t, 
using the time-dependent NEGF formalism, an adiabatic expansion in the slow time variable 
t = ^f- and Fourier transforming with respect to the fast variable [t — t') . The time- 

dependent nonequilibrium retarded and advanced Green's functions for the SWCN can be 
derived as Gc"(e, i) = [e - Hc(t) - Ss(e, t) - SD(e, where Ss(e, t) = hl^gs{e, i)hsc and 
SD(e, t)= /icDfl'D(e, ^)^cD self-energy terms due to metalhc contacts with g^s,D}i^, t) = 

|.-i/,,o,(£)l- the G.ee„-s toctio,. of the contacts. Th. . done by expanding the Gteen-. 
functions up to linear order in the slow variable t using the adiabatic expansion llll |: 

G'^^it -t',t) = -t',t}+ - t',t), where G(°>'"(t -t',t) = G'^^it - t',t)\i=t 

and G^'^^^'"'{t — t' ,t) = (^^) ^ (t — t', t) \t=t are the zeroth and first order Green's functions. 
Taking the Fourier transform with respect to the fast variable {t — t') the Green's functions 
become G'''"-{e,t) = G^^^'''''{e,t) + (^(^^'''"(e, t). We assume that the coupling to the contacts 
effectively gives rise to a finite imaginary term in the self-energies which is larger than 



IT] 



12| at all times. Therefore, we drop the term irj in the Green's function matrix for 



the SWCN. The coupling functions are considered to be energy independent Ts/D(t,t') = 
T{t,t') X exp[i /^*dtiAs/D(ti)] with r(t,t') = 271^,^^^ B Pa^oi,n{t)V*,^{t') |3|. These functions 
can also be expanded using the adiabatic expansion as r{s,D}(i — t',i) = ~ + 

rfl^.{t-t',t) withrf) (t-t',t)=r|s,D}(t-t',t)lf=t = rgD}(^l=2^E.es,DP«^",nW^;,nW 



as the zeroth and 



and rgj,^(t) = (f - t)A(t) = {f - t) [l'-^^\t=t - ir;^3^j3j(t)As/D(0) 
first order coupling functions. Using these nonequilibrium Green's functions and coupling 
functions we can identify the zeroth and first order contribution to the current: J{t) = 
J(0)(t) + J«(t). 
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We explicitly calculate the Green's function (Gc) for the Hamiltonian corresponding to 
:igure 2. In the matrix form, the Hamiltonian of the system can be divided into three blocks 



13| corresponding to the semiconducting SWCN and the two metallic contacts S and D 



m 



Hs hsc 

/iL Hn he 



\ 



hl^ ) 



(2) 



where -f^{s,D} are the contact Hamiltonians and Hq is the time-dependent tight-binding 
Hamiltonian for the SWCN with matrix elements: -ff^.^^=£:^j(t), HBiB^=£Bi{t) and 
HA,B=HB,A=lii{t) Q, HB,A,+i=HA,+^B=li,i+i{t), and the rest of the off diagonal 
elements are zero. The coupling matrices h^c and /icd will be non zero only for nearest- 
neighbors. Assuming the tight-binding model and using the above Hamiltonian the zeroth 
order time-dependent retarded Green's function G^^'^'^ can be written in an M x M matrix 
form as 



( e-EAA^-^^it) 711 (t) 

711© s-eB^) li2{t) 

712(f) e-EA-A^ 



(3) 



III. RESULTS AND DISCUSSION 



The zeroth order Green's functions (G'-^-'^, G*^'^-''', and G*^*^^^) lead to the zeroth order 
current J^^\ which is the Landauer formula , M with the slow time variable t. The 
Landauer formula is derived using the adiabatic expansion for the Green's functions and the 



coupling functions in the expression of the time-dependent current equation (6) of 
is found to be 



I and 



J''\i) = lj^ Tr{Tf{t)Gt^\e^i)T^^\^Gf\emh{e) - Ue)]. (4) 



The transmission function of the system is identified as T{e, t) = 
Tr{rP{t)G'^c^"{e,i)r'^\i)G'^c^^{e,t}), where GS?^"'"(£,t) are the zeroth order retarded 
and advanced Green's functions of the SWCN at i and /{s,d}(^) are the Fermi distribution 
functions in the source and drain contacts. The couphng functions are related to the 
self-energies by the relationship r{g_D}(f) = ip{s,D}(t) — Sjg D|(t)], where we have considered 
self-energies to be energy independent. S{s,d}(^) contains both real and imaginary parts. 

To calculate the transmission function, we have assumed that only the first ele- 
ment of the rg'^(f)(Ss(t)) matrix and the last element of the rQ^(t)(SD(t)) matrix 
are present. The rest of the elements of matrices are considered to be zero. Then, 
the transmission function depends only on one off diagonal element of Gc'. T{e,t) = 
-'^3,11 (^^cim(^; t). 



To get an explicit expression for the current, we consider linear response [12|, |13|. In linear 
response, equation (4) becomes J'^°^(t) = | / ^^'(e, f)(5[/(£:— /is)~/(^~/^D)] where /i{s,D} are 
the chemical potentials associated with the source and drain contacts. This equation then 
can be written as J''°\i) = |r(£/, f)[/xs - /xd] as 6[f{e-fj.s) - f{6- fir,)] = (/^s - AiD)(-|f ) 



and (-f ) = 6{ej 



In the experiment pi, there are two DNA sequences: sequence 1 with 21 bases and 
sequence 2 with 24 bases. Hence, we have two matrices, one is a 23 x 23 matrix and the 



2|. 



other is a 26 X 26 matrix 



16|. For the above Hamiltonian matrix and using rg'^|^(r) = 
-2Jm(Ss,ii(t)) = -2/mSs(F) and r|^^^j^j^(t) = -2/m(SD,MM(i)) = -2/mSD(t), the general 
expression for the transmission function of a SWCN decorated with DNA sequence 1 (2) 
and gas is given by 

where p = q = M/2 for even M and p = (M - l)/2 and q = (M + l)/2 for odd mQ. 
Equation (5) gives an explicit formula for the transmission function in terms of and 
indicating the dependence of the transmission function and hence the current on the hopping 
integrals and the on-site energies which are functions of time t. 

We also derive the first order contribution to the time-dependent current using equation 
(6) of Isl: 
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— —ImTr 
n 



|/s/D(.){r(°)^/°(t)G«^-iA(t)^ 



(6) 



where we have used t) = d(t - t')e'^(*-*')(t - t')G(°)'^'^(t - t',t) and the 

dependence of and G(^)'''" on e and i has been suppressed. An expression for 

G^^^''''"-(e, t) can be exphcitly calculated from G^'"'{t — t', t) and equation (3) using the formula: 

To analyze the sensor response in terms of the hopping integrals and on-site en- 
ergy, let us consider the model parameters jij and to have the form '~fij{Aacc{t)) = 
7oexp(=FAacc(t)/ao) and exX^'^ccii)) = £^oexp(=FAacc(t)/ao) where — (+) sign appears when 
ttcc decreases (increases). Here 70 and Eq are the hopping integral and on-site energy without 
the gas and DNA with ao = 0.33^4°. 7jj(Aacc(t)) and exX^^^ccii)) are the modified param- 
eters when the nearest-neighbor carbon-carbon distance Occ changes by Aocc as a function 
of time i due to the interaction of gas molecules with the DNA-decorated SWCN at each 
i s. As the experimental observations [5] show that the sensor response to PA differs in 
sign and magnitude from the response to other gases, the decreased sensor response [AI/Iq, 
table 1 ^) to all the gases except PA is due to —Aacdt) which enhances 7jj(Aacc(t)) and 
hence the electron current resulting in a decrease in the hole current (here we have assumed 
only the contribution of jij and neglect eXi in (5)). The increased sensor response to PA 
is because of +Aacc{t) which reduces 7jj(Aacc(t)) and the electron current and therefore 
increases the hole current. 

To demonstrate that the theory and experiment are in agreement we explicitly calculate 
the values for AJ^^^{t)/jQ (AI/Iq of [5]) for a 7 x 7 matrix (for simplicity of the calculation) 
using equations (4), (5) and the form for jij given in the above paragraph, for methanol with 
DNA sequence 2. Here Jo is the current without the odour. The table 1 given below gives 
the values for different parameters considered in the calculation. Figure 3 is the plot between 
the sensor response, calculated using the parameters given in table 1, and the time. This 
shows a good match between the experimental and theoretical results. Hence, the formula 
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reproduces the current characteristics of the experiment The results can be used to 
calculate the two time current correlation function and compared with the plots calculated 
from the experiment js], figure 1 and will be reported elsewhere. 

Using J{t) the charge transfer per second can be calculated as AQ{ti) = Jiti) — J(tj_i). 
These are some predictions of the model for the experiment. Using the expressions for G'-^^ 
and we find that J^^\t) changes with ^ij{t), exXt), ^sAt), and 



at 



18| while J'^°^(t) changes with jijii), exXi) and Ss,d(^- 



IV. CONCLUSIONS 

In conclusion, the two time current correlation function for the experimental data 5| has 
been calculated and for this system, analytical calculations of quantum electronic transport 
have been presented by setting up a tight-binding model and applying the time- dependent 
NEGF formalism. The Green's functions and coupling functions have been expanded using 
the adiabatic expansion in the slow variable and the Fourier transform has been taken with 
respect to the fast variable for the Green's functions. With the help of the Green's functions 
and coupling functions the zeroth and first order contributions to the current have been 
investigated. 

Equations (1), (4)-(6) carry the principle results of this paper. The correlation function 
equation (1) is a measure of the sensitivity and selectivity of the DNA-decorated SWCN 
sensors and suggest that these gas flow sensors may also be used as DNA sequence detectors 
where the pattern of correlation functions may be used as a benchmark for the particular 
chemical signal encoded in a DNA sequence. Equation (4) presents the zeroth order time- 
dependent current, which is the Landauer formula that depends on the slow time variable 
t. The dependence of the transmission function and hence the current on 7ij(t) and exXt) 
has been shown by an explicit formula, equation (5). Equation (6) shows the first order 
contribution to the current for the experiment jsl, which is proportional to the zeroth and 
first order Green's functions. The formula for the time- dependent current is then used to 
compare the theoretical results with the experiment. An expression for the charge transfer 
per second is obtained using the current. 

The numerical and analytical approaches used in this work can be applied to a broad 
range of systems where gas flows over nano-structures doped with different chemical and 



9 



biological molecules. This will provide a method for the study of time-dependent electronic 
transport in low dimensional disordered systems with gas flow. 
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t(s) t(s) 

FIG. 1: The two time current correlation function C{t) versus t(s) for (a) DMMP + sequence 2 
and DNT + sequence 1 (b) TMA + sequence 2 and PA + sequence 1 (c) methanol + sequence 1 
and (d) methanol + sequence 2. 



TABLE I: Values for different parameters used in the calculation of the current with 70=2. 5eV for 
methanol with DNA sequence 2. The units of ej and A 

dec 

eV and A respectively. 



tis) 


h 


t2 


^3 


t4 




Aj(°)(t)/Jo 


-0.0328 


-0.0703 


-0.0824 


-0.0776 


-0.1093 


Efit) 


0.01 


0.1 


1.4 


1.6 


2.8 




-0.004 


-0.004 


-0.004 


-0.004 


-0.004 


AaBiA2{i) 


-0.006 


-0.012 


-0.012 


-0.012 


-0.012 






-0.006 


-0.012 


-0.012 


-0.012 








-0.004 


-0.008 


-0.008 


AaA^Bait) 








-0.006 


-0.012 


AttBsAAt) 










-0.006 
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FIG. 2: Portion of an unwrapped SWCN. ^ijit) are the hopping integrals indicating hopping 
between different carbon atoms Ai and Bi at time i. Here, t = ti, • ' ' j ^Af- 7o represents the 
hopping integral in the absence of the gas and DNA. /isi represents the coupling between the source 
contact and the first carbon atom Ai while /ia/d is the coupling between the last carbon atom Xm 
and the drain contact. The ovals represent the DNA bases cytosine, thymine etc. while the circles 
indicate the gas molecules. 
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t(s) 

FIG. 3: The sensor response AJ^^\t)/ Jq versus time t(s) for methanol with DNA sequence 2. 
This plot shows agreement between the theory and experiment. Here, at ti the theoretical and 
experimental values of the current are different which is due to experimental error as this is the 
time when the gas molecule starts replacing the air molecule. We have found similar results for 
methanol with the other DNA sequence 1. Experimental data is used from with permission 
from American Chemical Society (2010). 



14 



